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Abstract 

The colour fields created by the static tetraquark and pentaquark systems are computed in 
quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24^ x48 lattice at ^5 = 6.2. 
We generate our quenched configurations with GPUs, and detail the respective benchmanrks in 
different SU(N) groups. While at smaller distances the coulomb potential is expected to dominate, 
at larger distances it is expected that fundamental flux tubes, similar to the flux-tube between 
a quark and an antiquark, emerge and confine the quarks. In order to minimize the potential 
the fundamental flux tubes should connect at 120o angles. We compute the square of the colour 
fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and with 
APE smearing. The tetraquark system is well described by a double- Y-shaped flux-tube, with two 
Steiner points, but when quark- antiquark pairs are close enough the two junctions collapse and we 
have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well described by a 
three- Y-shaped flux-tube where the three flux the junctions are Steiner points. 

1 Motivation 

Multiquark exotic hadrons like the tetraquark and the pentaquark, different from the the ordinary 
mesons and baryons, have been studied and searched for many years. The tetraquark was initially 
proposed by Jaffe [Ij as a bound state formed by two quarks and two antiquarks. Presently several 
observed resonances are tetraquark candidates. The most recent tetraquark candidates have been 
reported by the Belle Collaboration in May, the charged bottomonium (10610) and (10650) [2j. 
However a better understanding of tetraquarks is necessary to confirm or disprove the X, Y and possibly 
also light resonances candidates as tetraquark states. 

On the theoretical side, the first efforts have been to search for bound states below the strong decay 
threshold [21 HJ O |6] , as it is apparent that the absence of a potential barrier may produce a large decay 
width to any open channel. Recent investigations found that the presence of an angular momentum 
centrifugal barrier may increase the stability of the system [71 18] . 

In the last years, the static tetraquark potential has been studied in Lattice QCD computations 
[9l [ini [H]. The authors concluded that when the quark-quark are well separated from the antiquark- 
antiquark, the tetraquark potential is consistent with One Gluon Exchange Coulomb potentials plus a 
four-body confining potential, suggesting the formation of a double-Y flux tube, as in Fig. [T} composed 
of five linear fundamental flux tubes meeting in two Format points [HJ [12l |T3j. A Format, or Steiner, 
point is defined as a junction minimizing the total length of strings, where linear individual strings join 
at 120° angles. When a quark approaches an antiquark, the minimum potential changes to a sum of 
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Figure 1: We compare (a) the tetraquark flux tube (or string) model, the elementary flux tubes meet in 
two Fermat points, at an angle of a = 120° to form a double- Y flux tube (except when this is impossible 
and the flux tube is X-shaped) with (b) Lagrangian density 3D plot for ri = 8, r2 = 14, presented in 
lattice spacing units (colour online). 

two quark- antiquark potentials, which indicates a two meson state. This is consistent with the triple 
flip-flop potential, minimizing the length, with either tetraquark flux tubes or meson-meson flux tubes, 
of thin flux tubes connecting the different quarks or antiquarks [6| [8] . 

Here we study the colour flelds for the static tetraquark system [H], with the aim of observing the 
tetraquark flux tubes suggested by these static potential computations. The study of the colour flelds 
in a tetraquark is important to discriminate between different multi-quark Hamiltonian models. Unlike 
the colour flelds of simpler few-body systems, say mesons, baryons and hybrids, |T5l [IHl [I8], the 
tetraquark flelds have not been previously studied in lattice QCD. 



2 Computing Fields with the Wilson loop and the Plaquette 



To impose a static tetraquark, we utilize the respective Wilson loop [9l [lO] of Fig. [2| given by W4Q 
|Tr(Mii?i2M2Li2), where 
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The chromoelectric and chromomagnetic flelds on the lattice are given by the Wilson loop and 
plaquette expectation values, 
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where the jk indices of the plaquette complement the index i of the magnetic fleld, and where the 
plaquette at position r = (x, ^, z) is computed at t = T/2, 



P,. (r) = 1 - ^ReTr [u,{v)U,{v + ^)f/t(r + v)Ul{v) 



(3) 
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Figure 2: Tetraquark Wilson loop as defined by Alexandrou et al , and by Okiharu et al 

The energy (Ti) and lagrangian (C) densities are then computed from the fields, 

(^(r)> = l{{E\r)) + {B'ir))) , (4) 
(£(r)) = I ((E2(r)) - {B'{r))) . (5) 

To produce the results presented in this work , we use 1121 quenched configurations in a 24^ x 48 
lattice at /3 = 6.2. We also test that these configurations are already close to the continuum limit in a 
larger, 32^ x 64 lattice. We present our results in lattice spacing units of a, with a = 0.07261(85) fm or 
a-^ = 2718 ± 32 MeV. We generate our configurations in NVIDIA CPUs of the FERMI series (480, 580 
and Tesla 2070) with a SU(3) CUDA code upgraded from our SU(2) combination of Cabibbo-Marinari 
pseudoheatbath and over-relaxation algorithm [191 ED] . 

To compute the static field expectation value, we plot the expectation value {Ef{r)) or {Bf{r)) as a 
function of the temporal extent T of the Wilson loop. In order to improve the signal to noise ratio of the 
Wilson loop, we use 50 iterations of APE Smearing with w = 0.2 (as in [ill) in the spatial directions and 
one iteration of hypercubic blocking (HYP) in the temporal direction. with ai = 0.75, 0^2 = 0.6 
and 0^3 = 0.3. At sufficiently large T, the groundstate corresponding to the studied quantum numbers 
dominates, and the expectation value of the fields tends to a horizontal plateau, for each point r 
determined by the plaquette position. For the distances ri and r2 considered, we find in the range of 
T G [3, 12] in lattice units, horizontal plateaux with a /dof G [0.3, 2.0]. We finally compute the error 
bars of the fields with the jackknife method. 

3 The Tetraquark fields 

In our simulations, the quarks are fixed at (± ri/2, — r2/2, 0) and the antiquarks at (± ri/2, r2/2, 0), 
with ri extending up to 8 lattice spacing units and r2 extended up to 14 lattice spacing units, in order 
to include the relevant cases where r2 > V^ri. Notice that in the string picture, at the line r2 = V^ri in 
our (ri, r2) parameter space, the transition between the double-Y, or butterfiy, tetraquark geometry in 
Fig. [Ta|to the meson-meson geometry should occur. The results are presented only for the xy plane since 
the quarks are in this plane and the results with z are less interesting for this study. The fiux tube 
fields can be seen in Fig. [lb} |3] and |4} Theses figures exhibit clearly tetraquark double-Y, or butterfiy, 
shaped fiux tubes. The fiux tubes have a finite width, and are not infinitely thin as in the string models 
inspiring the Fermat points and the triple fiip-fiop potential, but nevertheless the junctions are close 
to the Fermat points, thus justifying the use of string models for the quark confinement in constituent 
quark models. 
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Figure 3: Lagrangian density for r2 = 14 and ri from to 6. The black dot points correspond to the 
Fermat points. The results are presented in lattice spacing units (colour online). 



In Fig. 5 (a) , we plot the chromoelectric field along the central fiux tube, \Ey^ at x = 0, for 
ri = 8, r2 = 14. As expected, the chromoelectric field along y is in agreement with the position of the 
Fermat points. The chromoelectric field along the x = central axis is maximal close to the Fermat 
points situated dit x —4.69 and ai x 4.69, fiattens in the middle of the fiux tube. Outside the fiux 
tube, the chromoelectric field is almost residual. 

In Fig. [5](b), we compare the chromoelectric field for the tetraquark and the quark-antiquark system 
in the middle of the fiux tube between the (di) quark and the (di)antiquark. As can be seen, for our 
larger distance r2 = 14 where the source effects are small, the chromoelectric field is identical up to the 
error bars, and this confirms that the tetraquark fiux tube is composed of a set of fundamental fiux 
tubes with Fermat junctions. 

To check which of the colour structures, tetraquark or meson-meson, produces the groundstate fiux 
tube, we study the x^/dof of the T plateaux. Surprisingly, event at distances as small as r2 — |riv^, 
where the fiip-fiop potential favours the two-meson fiux tube, we find T plateaux with a good /dof. 
This shows that the mixing between the tetraquark fiux tube and the meson-meson fiux tube is small, 
and it is possible to study clear tetraquark fiux tubes even at relatively small quark-antiquark distances. 
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Figure 4: Colour fields, energy density and Lagrangian density for ri = 8 and r2 = 14. The black 
dot points correspond to the Fermat points. The results are presented in lattice spacing units (colour 
online) . 

4 Foreword 

• The fiux tubes remain interesting in Lattice QCD, our results support the string model of con- 
finement, in particular for the tetraquark static potential 

• The mixing between the tetraquark and meson-meson fiux-tubes is small, which may contribute 
for narrower tetraquark resonances. 

• We are now studying in more detail other fiux tubes, as for the pentaquark in Fig. [6| 

• We will soon have codes for SU(2), SU(3), SU(4), etc, available on our webpage [121 ED]- 



TMs work was partly funded by the FCT contracts, PTDC/FIS/100968/2008, CERN/FP/109327/2009, 
CERN/FP/116333/2010 and CERN/FP/116383/2010, 
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Figure 5: (a) \Ey^ in the central axis x = for ri = 8, r2 = 14. We show with vertical dashed lines the 
location of the two Fermat points, (b) Profile cut at y = of the chromoelectric field for the tetraquark 
and quark-antiquark systems in the middle of the fiux tube. The results are presented in lattice spacing 
units (colour online). 
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